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Abstract 

A quantum mechanical search procedure to determine the real zeros of a poly- 
nomial is introduced. It is based on the construction of a spin observable whose 
eigenvalues coincide with the zeros of the polynomial. Subsequent quantum me- 
chanical measurements of the observable output directly the numerical values of 
the zeros. Performing the measurements is the only computational resource in- 
volved. 

PACS: 03.67.-a, 03.65Sq 
Introduction 

Quantum mechanical measurements are a computational resource. Various quantum 
algorithms use projective measurements at some stage or other to determine the period 
of a function e.g. [T]. In [2], projective measurements have been assigned a crucial role 
for a particular scheme of universal quantum computation which requires measurements 
on up to four qubits. Related schemes have been formulated based on measuring on 
triples and pairs [Sj, and finally on pairs of qubits only [I]. Measurements are also an 
essential part of Grover's search algorithm in order to actually read the result of the 
computation E] • 

A conceptually different strategy has been applied to propose a special-purpose ma- 
chine which is capable to diagonalize any finite-dimensional hermitean matrix by genuine 
quantum means, i.e. quantum measurements [3 IS]- In this approach of quantum diago- 
nalization, a hermitean matrix is considered as a quantum mechanical observable of an 
appropriate one-spin system. Projective measurements with a generalized Stern-Gerlach 
apparatus provide directly the unknown eigenvalues of the matrix which solves the hard 
part of the diagonalization. In this paper, it will be shown how to find the real zeros of 
a prescribed polynomial in a similar way, using quantum mechanical measurements. 



The quantum search procedure 



Consider a polynomial of degree N which is assumed to have N real zeros £, 



N 



P( x ) = ^2PnX n , p n e\R, PJV = 1- 



(1) 



n=0 



The assumption p^ = 1 is not a restriction since two polynomials Q(x) and Q(x)/c, c ^ 
0, have the same zeros. The quantum procedure to identify the zeros ( n of P(x) consists 
of two steps. First, one needs to find a hermitean companion matrix C of the poly- 
nomial P(x). By construction, its eigenvalues coincide with the unknown zeros of the 
polynomial P(x). Second, one determines the eigenvalues of the matrix C by a method 
inspired by the quantum diagonalization of a hermitean matrix. Effectively, they are 
obtained by measuring the eigenvalues of a quantum mechanical spin observable C with 
matrix representation C. 

Hermitean companion matrix of a polynomial 

It is straightforward to calculate the characteristic polynomial -Pm(A) of an (N x N) 
matrix M: 



where E is the (N x N) unit matrix. The polynomial -Pm(A) has degree N, and its zeros 
coincide with the eigenvalues /i n ,n = 1, . . . , N, of the matrix M. Hermitean matrices 
have real eigenvalues only, hence all zeros of its characteristic polynomial are located on 
the real axis. The inverse problem reads: 

• Given a polynomial P(x) of degree iV with real zeros, find a hermitean matrix C 
such that its characteristic polynomial is P{x). 

The matrix C is known as companion matrix of the polynomial P(x). Obviously, the 
companion matrix should be determined without reference to the roots of the given 
polynomial. If C is one solution of the inverse problem, then the matrices Cy = UCU^ 
provide solutions as well, where U is any unitary (N x N) matrix. The difficult part of 
the inverse problem lies in the requirement to find a hermitean matrix C: it is easy to 
specify a non-hermitean companion matrix Cq of the polynomial (0), viz., 



det(M-AE) =P M (A), 



(2) 
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Partial solutions of the problem to find hermitean companion matrices have been ob- 
tained in 0. A complete and constructive solution can be found in JU] 5 where a tridi- 
agonal companion matrix C is specified in terms of the coefficients p n , n = 0, 1, . . . , N. 
Explicitly, a polynomial P(x) of degree N with N real zeros can be written in the form 



P(x) = (-) N det (C - xE) , 
where the matrix C is tridiagonal and real symmetric, hence hermitean: 
/ 
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The nonnegative numbers dj.,k = 1, 2, . . . , N — 1, and the polynomials q n (x),n = 
1,2, ...,N, are generated when applying a Modified Euclidean Algorithm [Ej to the 
polynomial P(x). The operations required to determine the matrix elements of C are (1) 
repeated division of polynomials, (w) evaluation of specific polynomials at x = 0; (in) 
taking square roots of numbers c4- As a corollary, the Modified Euclidean Algorithm 
checks whether the given polynomial has real zeros only: if any of the numbers dk is 
found to be negative, P(x) can not have real zeros only. For completeness, the algorithm 
is sketched in the Appendix. 



Quantum search for eigenvalues of hermitean matrices 

Four steps are necessary to find the eigenvalues of a given hermitean matrix C with 
iV different eigenvalues by means of quantum measurements. Here an outline of this 
approach will be given only; for details about the procedure for (N x N) matrices readers 
should consult [7 , while it is illustrated for (2 x 2) matrices in [Hj. 

The matrix must be (2) written in standard form; next, it is interpreted as {11) 
matrix representation of a unique observable C of a quantum spin; this observable can be 
measured by (m) a specific apparatus which needs to be identified and built; finally, the 
apparatus is used to generate the eigenvalues by (iv) actually measuring the observable 
C. 

(1) Standard form of C: Write the hermitean (N x N) matrix C as a combination of 
linearly independent hermitean multipole operators J u , v — 0, . . . , iV 2 — 1, 

N 2 -l 

C= £ cT„, c„ = -Tr [CT„]elR. (6) 

u=0 iV 
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There are N 2 self-adjoint multipole operators T v = with matrix representa- 
tions T u , and they form a basis in the space of hermitean operators acting on an 
iV-dimensional Hilbert space 7i [T2|. These operators consist of all traceless sym- 
metric products of up to N spin components S = (S x , S y , S y ), plus the identity 
operator. Two multipoles are orthogonal with respect to a scalar product defined 



as the trace of their product: (1/N) Tr TJT v i 



6, 



(n) Identification of an observable: On the basis of the expansion © interpret the 
matrix C as representing an observable C for a spin with quantum number s = 
(N-l)/2: 

N 2 -l 

C = C(S) = J2 c » f ^ (7) 

thinking of the multipoles as functions of the spin components, T v = T U (S). 

{in) Setting up a measuring device for C : Swift and Wright have shown in ^2] how to 
devise a generalized Stern-Gerlach apparatus which measures any spin observable 
C. The construction generalizes the traditional Stern-Gerlach apparatus which 
measures the spin component e n ■ S along a direction specified by a unit vector 
e n . Setting up this device requires to create arbitrary static electric and magnetic 
fields in the laboratory, consistent with Maxwell's equations. The procedure is 
made explicit in Once constructed, the apparatus will split an incoming 

beam of particles with spin s into N = (2s + 1) subbeams corresponding to the 
eigenvalues of C. The working principle is equivalent to that of a standard Stern- 
Gerlach apparatus where C = S z . 

(iv) Determination of the eigenvalues: Prepare a spin s in a homogeneous mixture 
Po = I /(2s + 1). When carrying out measurements with the apparatus associated 
with C, the output of each individual measurement will be one of the eigenvalues 
( n of the matrix C. The actual values of the eigenvalues can be determined from 
the amount by which the particles in each subbeam are deflected from the straight 
line of flight. After sufficiently many repetitions, all eigenvalues ( n will be known. 
Since each eigenvalue occurs with probability 1/N, the probability not to obtain 
one of the N values decreases exponentially with the number of runs. 

By construction, the numbers ( n obtained in the last step coincide with the zeros of the 
polynomial P(x), and one can write 

N 

P(X) = Y[{X-Cn). (8) 

n=l 

The zeros ( n have been obtained by a genuinely quantum mechanical method. 
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Conclusions 

It has been shown that one can build a quantum mechanical special-purpose machine 
which is capable to output the roots of a polynomial of degree N. The underlying 
working principle is to perform appropriate quantum mechanical measurements using a 
generalized Stern-Gerlach apparatus. 

The possibility to extract information from individual quantum mechanical mea- 
surements may have implications for the interpretation of quantum mechanical states. 
Quantum root extraction of polynomials seems to strengthen the individual interpre- 
tation a single run of a measurement (with only one individual quantum system 
involved) provides information about one (randomly selected) zero of the polynomial; 
hence, no ensemble of identically prepared systems is required to obtain a useful answer 
from the experiment. 

Instead of repeating the experiment M ^> N times, one can imagine to run M 
identical experiments simultaneously, each generalized Stern-Gerlach apparatus being 
tuned to search for the zeros of the same polynomial P(x). If the input state is an 
M-fold direct product of the homogeneous mixture po, 

po ® • • • ® Po , (9) 

the resulting 'parallel' quantum search would, with large probability, produce all zeros 
( n at one go. In a sense, quantum mechanics is able to point almost instantaneously at 
the zeros of a given polynomial without any software program running. 

Appendix: Modified Euclidean Algorithm 

Given a polynomial 

P(X) =X JV +p JV _iX JV - 1 + ...+p , Pn £ I R , (10) 

the Modified Euclidean Algorithm defines recursively a number of polynomials 
Pk(x), k = 2, 3, . . . , iV — 1, of smaller degrees, and it generates other polynomials qk(x) 
and numbers dk which are required to define the hermitean companion matrix C. Start 
with 

P 1 (x)=P(x) P 2 ( X ) = , (11) 

and iterate the following steps. Divide Pk(x) by Pk+i(x), 

Pk{x) = q k {x)P k+ i(x) - r k (x) , (12) 

with a remainder r k (x) which may either be different from or equal to zero. Denote 
by d(Q(x)) the coefficient of the highest power of the polynomial Q(x): for example, 
d(P(x)) = 1. Define a polynomial Pk+2(x) and the number dk according to 
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W r k (x) ^ : P k+2 (x) = and d k = d(r k (x)) ; 

(m) r fc (x) = : P k+2 (x) = d ^^ x) and 4 = 0. 

The algorithm terminates if P k +i = 1, defining q k (x) = P k (x). Otherwise the proce- 
dure is repeated with P k+ i(x) and P k+ 2(x) in (fTTJ) and it will generate q k+ i and d k+ \, 
etc. 
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